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We present the recent robust determination of the value of the Dark Matter density at the Sun's 
location (p ) with a technique that does not rely on a global mass-modeling of the Galaxy. The 
method is based on the local equation of centrifugal equilibrium and depends on local and quite 
well known quantities such as the angular Sun's velocity, the disk to dark contribution to the 
circular velocity at the Sun, and the thin stellar disk scale length. This determination is inde- 
pendent of the shape of the dark matter density profile, the knowledge of the rotation curve 
at any radius, and the very uncertain bulge/disk/dark-halo mass decomposition. The result is: 
Pq = 0.43(0.1 l)(0.10)GeV/cm 3 , where the quoted uncertainties are due to the uncertainty a) in 
the slope of the circular-velocity at the Sun location and b) in the ratio between this radius and the 
exponential length scale of the stellar disk. The devised technique is also able to take into account 
any future improvement in the data relevant for the estimate. 
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1. Introduction 

Galaxy rotation curves (e.g. Rubin et al. 1980; Bosma et al. 1981) have unveiled the presence 
of a dark "mass component" in spirals. They are pillars of the paradigm of massive dark halos, 
composed of a still undetected kind of matter surrounding the luminous part of galaxies. The 
kinematics of spirals shows universal systematics (Persic, Salucci 1996; Salucci et al. 2007), which 
seems to be at variance with the predictions emerging from simulations performed in the A cold 
dark matter (ACDM) scenario, (e.g. Navarro et al. 1996), the current cosmological paradigm of 
galaxy formation (e.g. Gentile et al. 2004) (However, see Maccio et al. (2012)). 

Dedicated searches of DM particle candidates have intensified in recent years: direct-detection 
experiments look for the scattering of DM particles inside the detectors, which clearly is propor- 
tional to the DM density in the Sun's region p = Pdm{Rq)- Indirect-detection experiments search- 
ing for the secondary particles (e.g. neutrinos) produced by DM annihilation at the center of the 
Sun or Earth depend on the DM density inside these objects which in turn is driven, via a cap- 
ture mechanism, by the same local DM density p . Thus, in both direct and indirect searches the 
knowledge of the the local density p is very important for a precise prediction of the searched 
signal or to obtain reliable bounds on the DM particle cross- section. 

A value of p = 0.3GeV/cm 3 has been routinely quoted, but its origin is not clear, neither 
supported by data (see the introduction of Salucci et al. (2010)). An example is the work by 
Caldwell, Ostriker (1981) that devised what can be considered as the standard method to determine 
the value of p from a variety of observations. Their resulting value O^S+Q^GeV/cm 3 is however 
uncertain and moreover plagued by very outdated kinematics. 

In general, it is quite simple to infer the distribution of dark matter in spiral galaxies. Spiral's 
kinematics, in fact, reliably traces the underlying gravitational potential (Persic, Salucci 1996; 
Salucci et al. 2007). Then, from co-added and/or individual RCs, we can build suitable mass 
models that include stellar and gaseous disks, along with a spherical bulge and a dark halo. More 
in detail, by carefully analyzing (high quality) circular velocity curves, with the help of relevant 
photometric and HI data, one can derive the halo density at different radii. The accuracy of the 
measurements and of the analysis is excellent and the results obtained are at the heart of the present 
debates on Galaxy formation (e.g. Gentile et al. 2004, 2005; DeBlok 2009; an innovative review 
on this issue is found at www . sissa . it/ap/dmg/dmaw_presentation . html) 

To measure p , instead, is far from simple, because the MW kinematics, unlike that of external 
galaxies, does not trace the gravitational potential straightforwardly. We do not directly measure 
the circular velocity of stars and gas but rather the terminal velocity Vj of the rotating HI disk, 
and this only inside the solar circle (e.g. McClure-Griffiths, Dichey 2007). This velocity is related 
to the circular velocity V(r), for r < R & , by means of V(r) = Vj(r) +V r/R®, where R & ~ 8 
kpc is the distance of the Sun from the Galaxy center and V Q , the value of the circular velocity at 
the Sun's position. Both quantities are known within an uncertainty of 5% - 10% (e.g. McMillan, 
Binney 2009) which triggers a similar uncertainty in the derived amplitude and slope of the circular 
velocity. As a result, (see also Sofue 2009), the MW circular velocity from 2 kpc to 8 kpc is known 
within a not negligible uncertainty: 



V(r) = (215 ±30) km/s 

d\ogV {r) / d\ogr = a(r) =0.0 ±0.15. 



(1.1) 
(1.2) 
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The variations in equation are due to a mix of observational errors in the kinematics, to 
uncertainties in the values of R & and Vq and to actual radial variations of V. The first two trigger 



also part of the possible range of the velocity slope (1.2). Data show that the radial variations of 
a(r) are, instead, quite small: da{r)/dr ~ 0±0.03/kpc ~ 0; in the disk region the MW RC can 
be approximated by a straight line, whose slope however is mildly uncertain. 

The outer (out to 60 kpc) MW "effective" circular velocity V(r) = (GM(r) /r) 1 / 2 is even more 
uncertain and it depends on the assumptions made on dynamical and structural properties of its 
estimators. It appears to decline with radius, with quite an uncertain slope <ilogV(r)/<ilogr = 
-0.20+[j;^ (R G <r < 60kpc) (Battaglia 2005; Xue et al. 2008; Brown et al. 2009). These 
uncertainties, combined with the intrinsic "flatness" of the RC in the region specified above (that 
complicates the mass modeling even in the case of a high-quality RC (Tonini, Salucci 2004)), 
make it very difficult to obtain a reliable bulge/disk/halo mass model and consequently an accurate 
estimate of p . 

A different approach was devised by Salucci et al. (2010); the idea is to resort to the equation 
of centrifugal equilibrium (see Fall Efstathiou 1980 for details) and to use recent results found in 
external galaxies (Salucci et al. 2007). Let us start with 

V 2 /r = a H +a D + a B , (1.3) 

where an, ci£>, and a B are the radial accelerations generated by the halo, stellar disk, and bulge. 
Taking first the (quite good) approximation of spherical DM halo, we have an oc r~ 2 / r p//(?") r 2 dr. 



A similar relation holds for the bulge. Therefore, by differentiating equation (|1.3|), we obtain the 
DM density at any radius in terms of the local angular velocity co(r) =V/r, the RC slope a(r), the 
disk-to-dynamical mass ratio j8 (r) (see later), and the bulge density Ps(r): 

p H {r) = (D(r) 2 [F tot (a(r))-F D (f3(r))}-p B (r). (1.4) 

with F tot and F D known functions. 

In external spirals, equation ( |0| ) is not useful for determining the DM density at different 
radius because 1) it collapses for r < Rd where F tot ~ Fd and where the bulge density is dominating, 
Pb ^ Ph. and 2) the radial variations of a(r) have non-negligible observational uncertainties; 3) 
the quantity CO is known with far less accuracy than V, the main kinematical observational quantity. 

Instead, when estimating p Q , i.e. the density of the MW DM halo at a specific radius (the 
Sun position), the above drawbacks disappear, infact: 1) we have that F tot (R & ) » Fd(R & ) because 



Rq > 3Rd, so equation (1.4) does not collapse and, as a bonus, the most uncertain term of its r.h.s. 
turns out to be very small 2) eo Q is very precisely measured; and 3) at the Sun's position the bulge 
density Pb(Rq) i s totally negligible, < P///50 (e.g. Sofue et al. 2009). 

The method is clearly simpler for a spherically symmetric DM halo, and if we assume an 
innnitesimally thin disk for the distribution of stars in the Galaxy. Nevertheless, in Salucci et al. 
(2010) the effects of a possible halo oblateness and disk thickness were dealt with. 

Our claim (Salucci et al. 2010) is that by means of this method we to obtain an independent 
and reliable determination of the local DM halo density and of its intrinsic uncertainty. 



2. A new local method for determining p 







We model the Galaxy as composed by a stellar exponential thin disk plus an unspecified spher- 
ical DM halo with density profile p#(r). For the present work we can neglect the thick stellar and 
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the HI disk because their surface density, between 2kpc and R Q , are 100 to 5 times smaller than 
the stellar surface density (Nakanishi, Sofue 2003). Similarly, we neglect the stellar bulge because, 
as mentioned above, its spatial density at R Q is virtually zero (e.g. Sofue et al. 2009). It is worth 
noticing that the standard method of galaxy modeling and its variants cannot take these very sim- 
plifying assumptions because the global modeling involves all these mass components in a crucial 
way. 

Let us rewrite the equation of centrifugal equilibrium by subtracting the disk component from 
the total acceleration. From its radial derivative we then find 



Ph(t) 



X a 



AnGr 1 dr 



V 2 (r) 



(2.1) 



where X q is a factor correcting the spherical Gauss law used above in case of oblateness q of the 
DM halo (Salucci et al. 2010). The X q correction is very small: X q ~ 1.00-1.05. 

The disk component can be reliably modeled as a Freeman stellar exponential thin disk of 
length scale R D = (2.5 ±0.2) kpc (Picaud, Robin 2004; M. Juric et al. 2008; Robin et al. 2008; 
Reyle 2009). The stellar surface density is then: E(r) = (M D /2nR 2 :i )e^^ RD . Also, the disk can be 
considered infinitesimally thin. In fact, its thickness zo is small, zo ~ 250pc (M. Juric et al. 2008) 
and moreover z.q <C Rd < R&- We can thus write av(r) = G ^° 3 r (IoKo — I\K\)X ZQ , where I n and K„ 
are the modified Bessel functions computed at r/lR^, andX Z0 ~ 0.95 accounts for the nonzero disk 
thickness. 



Since only the first derivative of the circular velocity V(r) enters in ( |2.1[ ) and in any case this 
function in the solar neighborhood is almost linear, we can write 



V(r)=V Q [l + a Q (r-R Q )/R, 



©J) 



where a Q = cc(R Q ) is the velocity slope at the Sun's radius. Then equation ( |2.l| ) becomes 



4kG 



V 2 (r) 



(l + 2a 



H{r/Rv)X Z{ 



(2.2) 



(2.3) 



with 2H(r/R D ) = (3I Ko-IiK 1 ) + (r/R D )(IiK -I Ki). Equation (gj) holds at any radius outside 
the bulge region and measures Ph(R q ) = p by subtracting the density of the stellar component 
from the one of the whole gravitating matter. 

The disk mass can be parametrized (Persic, Salucci 1990) by Md = j8 1.1 G~ 1 VqRq, with 
j8 = V£/V 2 \r q , i.e. the fraction of the disc contribution to the circular velocity at the Sun. Finally, 
by exploiting the fact that the quantity V/R\r b = ft) = (30.3 ±0.3) km/s/kpc is measured with very 
high accuracy and much better than V & and R & separately (McMillan, Binney 2009; Reid 2009), 
we obtain 



p = 1.2x 10 



-27 § 



CO 



cm 3 \km/skpc 



X a 



[\+2a & )-\.\(3f(r 0D )X Z( 



(2.4) 



where r QD = R Q /Rv and /(r 0D ) = r 3 QD H(r QD ). 

Let us to focus on the advantages of this technique : a) it does not require assuming a particular 
DM halo density profile or the dynamical status of some distant tracers of the gravitational field; b) 
it is independent of the poorly known values of V & and of the RC slope at different radii ce(r); c) 
it does not depend on the structural properties of the bulge, which in the mass modeling, leads to a 
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degeneracy with the values of the mass of stellar disk and of the DM halo, d) it only mildly depends 
on the ratio r D, as well as, on the disk mass parameter p\ Let us also note that the determination 
depends on the RC slope at the Sun a , but in a well specified way. 



To proceed further we discuss the uncertainties on the parameters appearing in equation ([2.4]). 
Our determination depends on the ratio r@o = Rq/Rd- For this we adopt the reference value and 
uncertainty r Q o = 3.4 ±0.5, as suggested by the values of Rp discussed above and by the average 
of values of R G obtained recently: R Q = 8.2±0.5kpc (Ghez et al. (2008); Gillessen et al. (2008); 
Bovy et al. (2009)). This leads to f(r Q D) — 0.42 ±0.20, whose uncertainty propagates only mildly 



into the determination of p , because the second term of the r.h.s. of equation (2.4) is only one 
third of the first. 

Present data constrain the slope of the circular velocity at the Sun to a central value of a = 
and within a fairly narrow range —0.075 < a < 0.075. The uncertainty of a is the main source 
of the uncertainty of the present determination of p , see for instance (Oiling, Merrineld 2001). 

In equation ( |2.4] ), j8 is the only quantity that is not observed and therefore intrinsically un- 
certain. We can, however, constrain it by computing the maximum value j8 M for which the disk 
contribution at 2.2/?/) (where it has its maximum) totally accounts for the circular velocity. With 
no assumption on the halo density profile one gets j8 M = 0.85, independently of V and R G (Per- 
sic, Salucci 1990). However, this is really an absolute maximal value and it corresponds, out 
to R & , to a solid body halo profile: Vf, x R a '' with o& = 1. Instead, mass modeling performed 
so far for the MW and for external galaxies tend to find a lower value: G5/,(3/?d) < 0.8, which 
yields j3 M = 0.77. We can also set a lower limit for the disk mass, i.e. j3 m : first, the microlens- 
ing optical depth to Baade's Window constrains the baryonic matter within the solar circle to 
be greater than 3.91O 1O M (McMillan, Binney 2009). Second, the MW disk B-band luminos- 
ity L B = 2 x 1O 1O L and the very reasonable value M D /L B = 2 again imply M D ~ 4 10 10 M e . All 
this leads to j8" ! = j8 M /1.3 ~ 0.65. We thus take p = 0.72±[J$ as reference range. 

Using the reference values and expanding around their central values, we find 



P 3 = o.4 3 5fX 



cm J 



1 + 2.9 a e - 0.64 J8 - 0.72 + 0.45 r eD - 3.4 



0.1 I — -0.25 ) +0.10 { -7-0.95 J ±0.07 f % 30.3 

\kpc / V / \km/skpc 



. (2.5) 



This equation estimates the DM density at the Sun's location in an analytic way, in terms of the 
involved observational quantities at their present status of knowledge. The equation is written in a 
form such that, for the present reference values of these quantities, the term in the square brackets 
on the r.h.s equals 1, and the central result is p = 0.43GeV/cm 3 . As such, the determination is 



ready to account for future improved measurements by simply inserting them in the r.h.s. of (2.5) 



The next step is to estimate the uncertainty in the present determination of p . From equa- 
tion (2.5) and the allowed range of values discussed above, we realize that the main sources of 
uncertainty are ce Q , j8 and r o. The other parameters give at most variations of 2-3%, and can be 
neglected in the following. 

Then, first, it is illustrative to consider CC@, j8 and r £ as independent quantities. We thus have: 

/ \ GeV 

p = I 0.43 ±0.094 (as)T 0.016 (j3) ±O.O96 (r0fl) J — j , (2.6) 
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where Ai x \ means that A is the total effect due to the possible span of the quantity x. 

At this point, we can go one step further, and assume that the MW is a typical spiral, and using 
recent results for the distribution of matter in external galaxies, namely that DM halos around 
spirals are self similar (Salucci et al. 2007) so that the fractional amount of stellar matter directly 
dictates the value of rotation curve slope a (Persic, Salucci 1990): 

J8 =0.72 -0.95 a . (2.7) 



Using this relation in equation ( |2.5| ) we find (neglecting the irrelevant q and zq terms) 
GeV 



p = 0.43 



cm 3 



1 + 3.5 a e +0.451 r 0Z) -3.4 ) +0.07 ( - — " 30.3 



km/skpc 

From the current known uncertainties, with the estimated range of a & , we finally arrive to 

\GeV 



(2.8) 



/ \ tiev 

p = (0.430 ± 0. 113,,.. ± 0.096 , „ J — f . (2.9) 



Its uncertainty mainly reflects our poor knowledge of the velocity slope a G and the uncertainty 
in the galactocentric Sun distance. 

3. Discussion and conclusion 

We have described here a local determination of p , which relies directly on the equation of 
centrifugal equilibrium, by estimating the difference between the 'total' needed density and that of 
the stellar component. 

The method leads to a very reliable kinematical local determination of p , avoiding model- 
dependent and dubious tasks, mandatory with the standard method, i.e., a) to assume a particular 
DM density profile and a specific dynamical status for the tracers of the gravitational potential, b) to 
deal with the non-negligible uncertainties of the global MW kinematics, c) to uniquely disentangle 
the flatfish RC into the different bulge/disk/halo components. 

This estimate of p also shows that any kinematical determination using the Galactic rotation 
ultimately depends at least on three local quantities, the slope of the circular velocity at the Sun, the 
fraction of its amplitude due to the DM, and the ratio between the Sun galactocentric distance and 
the disk scale-length, whose uncertainty unavoidably propagates in the result. Two of these three 
quantities can be related by noting that the MW is a typical Spiral and using the relations available 
for these kind of galaxies (Salucci et al. 2007), so that the final uncertainty can be slightly reduced. 

The resulting local DM density that we find is p = (0.43 ± 0.11 (a; ) ±0.10( r , o ))GeV/cm 3 . 
We remark that it is free from theoretical assumptions and can be easily updated by means of 



equation ( |2.5| ) as the relevant observational quantities will become better known. 

Let us comment on other determinations in the literature. By applying a global modeling 
method with a refined statistical analysis to a large set of observational data, Catena, Ullio (2009) 
claimed a very precise measure, p = (0.389 ± 0.025 )GeV/cm 3 , in agreement with our results 
but at variance with our estimated uncertainty (see Salucci et al. (2010) and also Weber, de Boer 
(2009)). 

Finally, it is worth discussing the claim by Bidin et al. (2012). They, by means of the complex 
and not unambiguously accepted technique of exploiting the vertical velocity dispersion of old 
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tracer stars of the thick disk at the solar neighborhood, claimed the absence of any DM at the 
solar position, in disagreement with our (and others) result. However, as pointed out by Bovy 
and Tremaine (2012) their result is plagued by arbitrary assumptions: different z dependence of 
the tracers circular velocity and/or different kinematical status and spatial distribution can trigger a 
large number of different possible mass models, some of them clear of DM, others with proportions 
larger than those we estimate in (Salucci et al. 2010). As a result, the Bidin et al. (2012) method 
needs a series of independent checks and increased precision/statistics before being able to deliver 
a robust determination of p . 

In conclusion, we believe that our technique provides the most trustable estimate of the local 
dark matter density and of its uncertainties. 
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